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ABSTRACT 

The Kepler Mission has recently announced the discovery of Kepler-lOb, the smallest exoplanet 
discovered to date and the first rocky planet found by the spacecraft. A second, 45-day period transit- 
like signal present in the photometry from the first eight months of data could not be confirmed as 
being caused by a planet at the time of that announcement. Here we apply the light-curve modeling 
technique known as BLENDER to explore the possibility that the signal might be due to an astrophysical 
false positive (blend) . To aid in this analysis we report the observation of two transits with the Spitzer 
Space Telescope at 4.5 /mi. When combined they yield a transit depth of 344±85 ppm that is consistent 
with the depth in the Kepler passband (376 ±9 ppm, ignoring limb darkening), which rules out blends 
with an eclipsing binary of a significantly different color than the target. Using these observations 
along with other constraints from high-resolution imaging and spectroscopy we are able to exclude 
the vast majority of possible false positives. We assess the likelihood of the remaining blends, and 
arrive conservatively at a false alarm rate ofl.6x ICC 5 that is small enough to validate the candidate 
as a planet (designated Kepler-lOc) with a very high level of confidence. The radius of this object 
is measured to be R p = 2.227^®'^ R§ (in which the error includes the uncertainty in the stellar 
properties), but currently available radial- velocity measurements only place an upper limit on its 
mass of about 20 M®. Kepler-lOc represents another example (with Kepler-9d and Kepler-llg) of 
statistical "validation" of a transiting exoplanet, as opposed to the usual "confirmation" that can take 
place when the Doppler signal is detected or transit timing variations are measured. It is anticipated 
that many of Kepler 's smaller candidates will receive a similar treatment since dynamical confirmation 
may be difficult or impractical with the sensitivity of current instrumentation. 

Subject headings: binaries: eclipsing — planetary systems — stars: individual (Kepler-10, KOI-072, 
KIC 11904151) — stars: statistics 



1. INTRODUCTION 
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The Kepler Mission has recently made public a cata- 
log of all transiting planet candidates identified during 
the first four months of observation by the spacecraft 
(jBorucki et all 1201 lbt l. Included in this list of 1235 ob- 
jects are nearly 300 in the category of super-Earths (de- 
fined here as having radii in the range 1.25 R® < R p < 
2Rq), and several dozen of Earth size (R p < 1.25 -R©). 
The wealth of new information promises to revolutionize 
our knowledge of extrasolar planets. Although strictly 
speaking these are still only candidates since confirma- 
tion by spectroscopic or other means is not yet in hand, 
expectations are high that the r ate of false positives 
in this list is relatively s mall (see iBorucki et al.l 1201 lbt 
iMorton fc Johnson] 1 2 1 it ) . Consequently, results from 
this sample concerning the general properties of exo- 
planets have already begun to emerge, including studies 
of the ar chitecture and dynam ics of multiple transiting 
systems (|Lissauer et al.lr2011bf) . an inv estigation of the 
statis tical distribution of eccentricities dMoorhca d et al.l 
1201 If ), and first estimates of the rate of occurrence of 
plane ts larger than 2i? m with orbital periods up to 50 
days (Howard et al.N20fll ). among others. 

For good reasons the confirmation or "validation" of 
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small transiting planets 18 (Earth-size or super-Earth- 
size) has attracted considerable attention, but has proven 
to be non-trivial in many cases because of the difficulty 
of detecting the tiny radial- velocity (RV) signatures that 
these objects cause on their parent stars , as ex emplified 
by t he cases of CoRoT -7b (jLeger et al J 12009ft . Kepler- 
9d dTorres et al.l 120 111) , and Kepler-llg (|Lissauer et al.l 
l2011al) . In fact, such spectroscopic signals are often too 
small to detect with current instrumentation, and the 
planetary nature of the candidate must be established 
statistically, as in the latter two cases. 

The smallest planet disc overed to date , Kep ler-10 b, 
was announced recently by iBatalha et al.l ((20111 ) , and is 
the Kepler Mission's first rocky planet. It has a mea- 
sured radius of 1.416lo 036 ^® anc ^ a mass 01 4.6+13 -^ffi> 
leading to a mean density of 8.8+2'gg cm -3 that im- 
plies a significant iron mass fraction (Batalha et al. 
I20P1 . Its parent star, Kepler-10 (KIC 11904151, 
2MASS 119024305+5014286), is relatively bright among 
the Kepler targets (Kepler magnitude Kp = 10.96) and 
displays two periodic signals with periods of 0.84 days 
and 45.3 days, and flux decrements (ignoring limb dark- 
ening) of 152 ± 4 pp m and 376 ± 9 ppm, respectively 
(jBatalha et al.|[201lD . The extensive observations that 
followed the detection of these signals are documented 
in detail by those authors, and include the difficult mea- 
surement of the reflex radial- velocity motion of the star 
with a semi- amplitude of only 3.3]t?Qiiis — 1 and a pe- 
riod that is consistent with the shorter signal. As is 
customary also in ground-based searches for transiting 
planets, the shapes of the spectral lines were examined 
carefully to rule out changes of similar amplitude cor- 
relating with orbital phase that might indicate a false 
positive, such as a background eclipsing binary (EB) 
blended with the target, or an EB physically associated 
with it. However the precision of the measurements (bi- 
sector spans) compared to the small RV amplitude did 
not allow such changes to be ruled out unambiguously. 
False positive scenarios were explored with the aid of 
BLENDER, a technique that models the transi t light curves 
to tes t a wide range of blend configurations (|Torres et al.l 
l20"H . and it was found that the overwhelming major- 
ity of them can be rejected. Th is and other evidence 
presented by IBatalha et al.l (|201lD allowed the planetary 
nature of Kepler-10 b to be established with very high 
confidence. 

This was not the case, however, for the 45-day period 
signal referred to as KOI-072.02 (Kepler Object of In- 
terest 72.02), which is the subject of this paper. No sig- 
nificant RV signal was detected at this period, and only 
an upper limit on its amplitud e could be placed. Using 
BLENDER, IBatalha et al.l (|201lD were able to rule out a 
large fraction of the blend scenarios involving circular or- 
bits (including hierarchical triples), but eccentric orbits 
were not explored because of the increased complexity 
of the problem and the much larger space of parameters 
for false positives. While circular orbits are a reasonable 

18 In the context of this paper "confirmation" refers to the unam- 
biguous detection of the gravitational influence of the planet on its 
host star or on other bodies in the system (e.g., the Doppler signal, 
or transit timing variations) to establish the planetary nature of 
the candidate; when this is not possible, we speak of "validation" , 
which involves an estimate of the false alarm probability. 



assumption for Kepler-10 b because of the strong effects 
of tidal forces at close range, this is not true for KOI- 
072. 02 on accoun t of its much longer orbital period (see, 
e.g., [Mazch 200j|); eccentric orbits can not be ruled out. 

This provides the motivation for the present work, in 
which we set out to examine all viable astrophysical false 
positive scenarios for KOI-072.02 with the goal of validat- 
ing it as a bona- fide planet. In addition to improvements 
in the BLENDER modeling, we bring to bear new near- 
infrared observations obtained with the Spitzer Space 
Telescope in which the transits are clearly detected, as 
well as the complete arsenal of follow-up observations 
gathered by the Kepler team, including high-resolution 
adaptive optics imaging and speckle interferometry, high- 
resolution spectroscopy, and an analysis based on the 
Kepler observations themselves of the difference images 
in and out of transit for positional displacements (cen- 
troid motion). All of these observations combined with 
the strong constraints provided by BLENDER significantly 
limit the kinds of blends that remain possible, and as 
we describe below they allow us to claim with very high 
confidence that KOI-072.02 is indeed a planet. Its esti- 
mated radius is approximately 60% of that of Neptune. 
With this, Kepler-10 becomes the Mission's third con- 
fir med multi-planet sys t em (after Kepler-9 an d Kepler- 
11: iHolman et al.l 120101: iLissauer et atl l2011af) contain- 
ing a transiting super-Earth-size planet and at least one 
larger planet that also transits. 

We begin with a brief recapitulation of the BLENDER 
technique, including recent improvements. We then 
present the Warm Spitzer observations at 4.5 /jm that 
help rule out many blends, and we summarize additional 
constraints available from other observations. This is 
followed by the application of BLENDER to KOI-072.02 in 
order to identify all blends scenarios that can mimic the 
Kepler transit light curve. Next we combine this infor- 
mation with the other constraints and carry out a sta- 
tistical assessment of the false alarm rate for the planet 
hypothesis, leading to the validation of the candidate as 
Kepler-10 c. We conclude with a discussion of the possi- 
ble constitution of the new planet in the light of current 
models, and the significance of this type of validation. 

2. REJECTING FALSE POSITIVES WITH BLENDER 

The detailed morphology of a transit light curve 
(length of ingress/egress, total duration) contains impor- 
tant information that can be used to reject many false 
positive scenarios producing brightness variations that 
do not quite have the right shape, even though they 
may well match th e observed transit depth (see, e.g. , 
iSnellen et all 12009). BLENDER (|Torres et al.ll200l (20111) 
takes advantage of this to explore a very large range of 
scenarios, including background or foreground eclipsing 
binaries blended with the target, as well as eclipsing bina- 
ries physically associated with the target in a hierarchical 
tri ple configura t ion. Following the notation introduced 
bv lTorres etaLl (|2011l) . the objects composing the binary 
are referred to as the "secondary" and "tertiary" , and the 
candidate is the "primary" . The tertiary can be either a 
star (including a white dwarf) or a planet, and the sec- 
ondary can be a main-sequence star or a (background) 
giant. 

With the help of model isochrones to set the stel- 
lar properties, BLENDER simulates blend light curves re- 
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suiting from the flux of the eclipsing pair diluted by 
the brighter target (and any additional stars that may 
fall within the photometric aperture). Each simulated 
light curve is compared with the Kepler observations in 
a x 2 sense to identify which of them result in accept- 
able fits (to be defined later). The parameters varied 
during the simulations are the mass of the secondary 
star (M 2 ), the mass of the tertiary (M 3 , or its radius 
i?3 if a planet), the impact parameter (b), the relative 
linear distance (d) between the eclipsing pair and the 
target, and the relative duration (D/D C i VC ) of the tran- 
sit compared to the duration for a circular orbit (see 
below). For convenience the relative linear distance is 
parametrized in terms of the difference in distance mod- 
ulus, AS, where AS = 51og(dEB/dKoi)- In the case of 
hierarchical triple configurations the isochrone for the bi- 
nary is assumed to be the same as for the primary (metal- 
lici ty of [Fe/H] = — 0.1 5 and a nominal age of 11.9 Gyr; 
see lBatalha et al1l201lD . whereas for background blends 
we have adopted for the binary a representative 3 Gyr 
isochrone of solar metallicity, although these parameters 
have a minimal impact on the results. For full details 
of the technique we refer the reader to the references 
above. Three recent changes and improvements that are 
especially relevant to the application to KOI-072.02 are 
described next: 

(i) The relatively long orbital period of KOI-072.02 
(45.3 days) precludes us from assuming that the eccen- 
tricity (e) is zero, as we were able to suppose in previous 
applications of BLENDER to Kepler-9d and Kepler-lOb, 
which have periods of 1.59 and 0.84 days, respectively. 
The reason this matters is that the duration of the tran- 
sit is set, among other factors, by the size of the sec- 
ondary star. Eccentricity can alter the speed of the ter- 
tiary around the secondary, making it slower or faster 
than in the circular case depending on the orientation 
of the orbit (longitude of periastron, uj). Given a fixed 
(measured) duration, blends with smaller or larger sec- 
ondary stars than in the circular case may still provide 
satisfactory fits to the light curve, effectively increasing 
the pool of potential false positives. 

BLENDER now takes this into account, although rather 
than using as parameters e and uj, which are the nat- 
ural variables employed in the binary light -curve gener- 
ating routine at the core of BLENDER (see iTorres et al.l 
a more convenient variable that captures the ef- 
fects of both i s the duratio n relative to a circular or- 
bit. Following iWinnl ([2010), this may be expressed as 
D/D c i IC fa yl — e 2 / (1 + esinw). Operationally, then, 
we vary D / D C1TC over wide ranges as we explore different 
blend scenarios, and for each value we infer the corre- 
sponding values of e and uj. In practice, in order to solve 
for {e, uj} from D/D c i IC it is only necessary to consider 
the limiting cases with to = 90° and 270°, corresponding 
to transits occurring at periastron and apastron, respec- 
tively, since these are the orientations resulting in the 
minimum and maximum durations for a given eccentric- 
ity. Other combinations of e and uj will lead to interme- 
diate relative durations that are already sampled in our 
D/D c i IC grid. It is worth noting that use of only these two 
values of uj leads to predicted secondary eclipses in the 
simulated light curves that are always located at phase 
0.5, whereas secondary eclipses in the real data might be 



present at any phase. For our purposes this is of no con- 
sequence, as KOI-072.02 has already had its light curve 
screened for secondary eclipses at any phase that might 
betray a false positive, as part of the vetting process. 
No such features are present down to the 100 ppm level. 
Thus, any simulated light curves from BLENDER that dis- 
play a significant secondary eclipse will yield poor fits no 
matter where the secondary eclipse happens to be, and 
will lead to the rejection of that particular blend scenario. 

(ii) For each false positive configuration BLENDER can 
predict the overall photometric color of the blend, for 
comparison with the measured color index of the can- 
didate as reported in the Kepler Input Catalog (KIC; 
iBrown et al.l I201TT) . A color index such as Kp — K s , 
where Kp is the Kepler magnitude and K s derives from 
the 2MASS catalog, provides a reasonable compromise 
between wavelength leverage and the precision of the in- 
dex. The latter varies typically between 0.015 and 0.030 
mag, depen ding on the passba nd and the brightness of 
the star (see IBrown et al.|[2011f) . We consider a particu- 
lar blend to be rejected when its predicted color deviates 
from the KIC value by more than three times the error of 
the latter. As it turns out, color is a particularly effective 
way of rejecting blends that include secondary stars of a 
different spectral type than the primary, such as those 
that become possible when allowing for eccentric orbits. 

(Hi) Recent refinements in the resolution of the 
BLENDER simulations to better explore parameter space, 
in addition to the inclusion of eccentricity (or D / D clrc ) as 
an extra variable, have increased the complexity of the 
problem as well as the computing time (by nearly two 
orders of magnitude) compared to the relatively simple 
case of circular orbits. The number of different parameter 
combinations examined with BLENDER (and correspond- 
ing light-curve fits) can approach 7 x 10 8 in some cases. 
Consequently the simulations are now performed on the 
Pleiades cluster at the NASA Advanced Supercomputing 
Division, located at the Ames Research Center (Califor- 
nia), typically on 1024 processors running in parallel. For 
convenience hierarchical triple configurations (4 parame- 
ters) and background/foreground blends (5 parameters) 
are studied separately, each for the two separate cases of 
stellar and planetary tertiaries (for a total of four grids) . 
One additional fit is carried out using a true transiting 
planet model to provide a reference for the quality of the 
false positive fits in the other grids. 

The discriminating value of the shape information con- 
tained in the light curves, mentioned at the beginning of 
this section, is highlighted b y our BLENDER result s for 
Kepler-lOb, as described bv iBatalha et al.l ()2011l ). In 
that study it was found that all background eclipsing bi- 
nary configurations with stellar tertiaries yield very poor 
fits to the Kepler light curve, and are easily rejected. 
The underlying reason is that all such blend models pre- 
dict obvious brightness changes out of eclipse (ellipsoidal 
variations) with an amplitude that is not seen in the 
data, and that are a consequence of the very short or- 
bital period. 19 Hierarchical triple scenarios were also 

19 Note that the present post-processing of Kepler data in prepa- 
ration for the BLENDER analyses (see Sect. [4{ artificially suppresses 
out-of-eclipse variations to some extent, typically by median filter- 
ing, so that the light curves for periods as short as that of Kepler- 
10 b (0.84 days) are rendered essentially flat except for the transits 
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excluded based on joint constraints from BLENDER and 
other follow-up observations. The only configurations 
providing suitable alternatives to the true planet sce- 
nario involved stars in the foreground or background of 
the target that are orbited by a larger transiting planet. 
The considerable reduction in the blend frequency from 
the exclusion of all background eclipsing binaries led to 
a false alarm probability low enough to validate Kepler- 
10 b with a very high level of confidence, independently 
of any spectroscopic evidence. This remarkable result 
speaks to the power of BLENDER when combined with all 
other observational constraints. It also assumes consid- 
erable significance for Kepler-lOb, given that it was not 
possible to provide s eparate proof of the p lanetary nature 
of this signal in the Batalha et al.l (|2011l ) study from an 
examination of the bisector spans. The scatter of the bi- 
sector span measurements (10.5 ms _1 ) was three times 
larger than the RV semi-amplitude (3.3 ms _1 ), rendering 
them inconclusive. 

The situation regardi ng the BLENDER anal ysis of the 
KOI-072.02 signal in the lBatalha et al.l (|2011[ ) study was 
very different: the orbital period is much longer, and 
ellipsoidal variations are predicted to be negligible, so 
that background eclipsing binaries with stellar tertiarics 
remain viable blends. This, and the added complication 
from eccentric orbits, hindered the efforts of those au- 
thors to validate this candidate. With the benefit of the 
enhancements in BLENDER described above, we are now 
in a better position to approach this problem anew. 

As follow-up observations provide important con- 
straints that are complementary to those supplied by 
BLENDER, and play an important role in determining the 
false alarm rate for the planetary nature of KOI-072.02 
(Sect. [S]), we describe those first below, beginning with 
our new near-infrared Spitzer observations. 

3. OBSERVATIONAL CONSTRAINTS 
3.1. Warm Spitzer observations of KOI-072.02 

KOI-072.02 was observed during two transits with 
the IRAC instrument on the Spitzer Space Telescope 
(|Werner et all l2004t iFazio et al.l 12004ft at 4.5 /zm (pro- 
gram ID 60028). The observations were obtained on UT 
2010 August 30 and November 15, with each visit last- 
ing approximately 15 hr lOmin. The data were gathered 
in full-frame mode (256 x 256 pixels) with an exposure 
time of 6.0 s per image, which resulted in approximately 
a 7.1s cadence and yielded 7700 images per visit. 

The method we used to produce photometric t i me se - 
ries from the images is described bv lDesert et al.l (2009). 
It consists of finding the centroid position of the stel- 
lar point spread function (PSF) and performing aper- 
ture photometry using a circular aperture on individual 
exposures. The images used are the Basic Calibrated 
Data (BCD) delivered by the Spitzer archive. These files 
are corrected for dark current, flat-fielding, and detector 

themselves. In this sense the situation is similar to that mentioned 
earlier regarding the presence of secondary eclipses: obvious ellip- 
soidal variability in the raw data would normally trigger a false 
positive warning during the vetting process, preventing the target 
from becoming an object of interest. But if it reaches KOI status, 
we assume that out-of-eclipse modulations are insignificant so that 
the comparison with any BLENDER model in which those variations 
are present is meaningful and would yield a poor fit, sufficient in 
most cases to reject the blend. 



non-linearity, and are converted to flux units. We con- 
verted the pixel intensities to electrons using the infor- 
mation on the detector gain and exposure time provided 
in the FITS headers. This facilitates the evaluation of 
the photometric errors. We extracted the UTC-based 
Julian date for each image from the FITS header (key- 
word DATE_OBS) and corrected to mid-exposure. We 
converted to TDB-based barycentric Jul ian dates using 
the UTC2BJD 20 procedure developed by lEastman et all 
ffi)Wi\ . This program uses the JPL Horizons ephemeris 
to estimate the position of the Spitzer spacecraft dur- 
ing the observations. We then corrected for transient 
pixels in each individual image using a 20-point sliding 
median filter of the pixel intensity versus time. To do 
so, we compared each pixel's intensity to the median of 
the 10 preceding and 10 following exposures at the same 
pixel position, and we replaced outliers greater than 4cr 
with their median value. The fraction of all pixels we 
corrected is 0.02% for the first visit and 0.06% for the 
second. 

The centroid position of the stellar PSF was deter- 
mined using the DAOPHOT- related procedures GCNTRD, 
from the IDL Astronomy Library 21 . We applied the APER 
routine to perform aperture photometry with a circular 
aperture of variable radius, using a range of radii between 
1.5 and 8 pixels in steps of 0.5. The propagated uncer- 
tainties were derived as a function of the aperture radius, 
and we adopted the aperture providing the smallest er- 
rors. We found that the transit depths and errors varied 
only weakly with aperture radius for all light-curves an- 
alyzed in this project. The optimal aperture was found 
to have a radius of 4.0 pixels. 

We estimated the background by examining a his- 
togram of counts from the full array. We fit a Gaussian 
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Fig. 1. — Spitzer transit light-curves of KOI-072.02 observed in 
the IRAC band-pass at 4.5 /an. Top: Raw measurements (black 
points) with the same data binned by two superimposed (12 s bins, 
red points). Bottom: Measurements combined from the two vis- 
its and binned in 36min bins (295 points per bin), along with 
the best-fit limb-darkened transit model (integrated over the same 
duration). Both the data and the model shown here have been 
corrected for instrumental errors. 



20 http: //astroutils . astronomy . ohio-state . edu/time/ 

21 http: //idlastro . gsfc .nasa.gov/homepage .html 
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curve to the central region of this distribution (ignoring 
bins with high counts, which correspond to pixels con- 
taining stars), and we adopted the center of this Gaus- 
sian as the value of the residual background intensity. 
As seen already in previous Warm Spit zer observations 
(jDeming et al.lboill iBeerer et all 1201 If ), we found that 
the background varies by 20% between three distinct lev- 
els from image to image, and displays a ramp-like behav- 
ior as function of time. The contribution of the back- 
ground to the total flux from the stars is low for both 
observations, from 0.1% to 0.55% depending on the im- 
age. Therefore, photometric errors are not dominated by 
fluctuations in the background. We used a sliding me- 
dian filter to select and trim outliers in flux and position 
greater than 5<r, representing 1.6% and 1.3% of the data 
for the first and second visits, respectively. We also dis- 
carded the first half-hour's worth of observations, which 
is affected by significant telescope jitter before stabiliza- 
tion. The final number of photometric measurements 
used is 7277 and 7362. 

The raw time series are presented in the top panel of 
Figure [TJ We find that the point-to-point scatter in the 
photometry gives a typical signal-to- noise ratio (S/N) of 
280 per image, which corresponds to 90% of the theoret- 
ical signal-to-noise. Therefore, the noise is dominated by 
Poisson statistics. 

3.2. Analysis of the Warm Spitzer light curves, and 

results 

In order to determine the transit parameters and as- 
sociated uncertainties from the Spitzer time series we 
used a transit light curve model multiplied by instrumen- 
tal deco rrelation functions, as described by IDesert et al.l 
(|2011al ) . The transit light curves were co mputed with the 
IDL t ransit routine OCCULTSMALL from iMandel fc Agoll 
(2002). For the present case we allowed for a single free 
parameter in the model, which is the planet-to-star ra- 
dius ratio R p /R+ (or equivalently, the depth, in the ab- 
sence of limb darkening). The normalized orbital semi- 
major axis (system scale) a/i?*, the impact parameter 
b, the period P, and the time of mid transit T c were 
held fixed at the val ues derived from the Kepler light 
curve, as reported by iBatalha et al.l (|2011l) and summa- 
rized below in Sect. [7] Limb darkening is small at 4.5 /im, 
but was nevertheles s included in our modeling using the 
4-parameter law bv lClaretl (120001) and theoretical coeffi- 
cients published by iSingl (|2010l ) . 

The Spitzer /IRAC photometry is known to be sys- 
tematical ly affected by the so-called "pixel- phase effect " 
(see, e.g.. lCharbonneau et aLll2005tlKnutson et al.ll2008l) . 
This effect is seen as oscillations in the measured fluxes 
with a period corresponding to that of the telescope 
pointing jitter. For the first visit this period was 70 min, 
and the amplitude of the oscillations was approximately 
2% peak-to-peak; for the second visit the period was 
35 min, and the amplitude about 1%. We decorrelated 
our signal in each channel using a linear function of time 
for the baseline (two parameters) and a quadratic func- 
tion of the PSF position (four parameters) to correct 
the data for each channel. We performed a simultane- 
ous Levenberg-Ma rquardt least-squares fit to the data 
(|Markward t 2009) to determine the transit and instru- 
mental model parameters (7 in total). The errors on each 
photometric point were assumed to be identical, and were 



set to the rms residual of the initial best fit. To obtain 
an estimate of the correlated and s ystematic errors in 
our measurements (|Pont et al.l l2006f we used the resid- 
ual permutation bootstra p technique, o r "Pra yer Bead" 
method, as described by IDesert et al.l {2009). In this 
method the residuals of the initial fit are shifted system- 
atically and sequentially by one frame, and then added 
to the transit light curve model before fitting again. We 
considered asymmetric error bars spanning 34% of the 
points above and below the median of the distributions 
to derive t he la uncertainties fo r each parameter, as de- 
scribed by IDesert et all (|2011bQ . 

The bottom panel of Figure [T] shows the best-fit model 
superimposed on the observations from the two visits 
combined, with the data binned in 36 min bins for clar- 
ity (295 points per bin). The transit depths at 4.5 /^m 
(after removing limb- darkening effects) are 353^3 ppm 

for the first visit and 339j^ for the second, which are 
in good agreement with each other. The weighted aver- 
age depth of 344 ± 85 is consistent with the non-limb- 
darkened va lue of 376 ± 9 ppm derived from the Kepler 
light curve (|Batalha et al.ll2011f) well within the lcr er- 
rors, strongly suggesting the transit is achromatic, as ex- 
pected for a planet. 

The above Spitzer observations provide a useful con- 
straint on the kinds of false positives (blends) that may 
be mimicking the KOI-072.02 signal. For example, if 
Kepler- 10 were blended with a faint unresolved back- 
ground eclipsing binary of much later spectral type that 
manages to reproduce the transit depth in the Kepler 
passband, the predicted depth at 4.5 fim may be expected 
to be larger because of the higher flux of the contaminat- 
ing binary at longer wavelengths compared to Kepler- 10. 
Since the transit depth we measure in the near infrared 
is about the same as in the optical, this argues against 
blends composed of stars of much later spectral type. 
Based on model isochrones and the properties of the tar- 
get star (see below), we determine an upper limit to the 
secondary masses of 0.77 Mq. This Spitzer constraint is 
used in Sect. 2] to eliminate many blends. 

3.3. Additional observational constraints on possible 
false positives 

Further constraints of a different kind are provided by 
high- resolution imag ing as described in more detail by 
IBatalha et al.l (| 2 1 If ) . Briefly, these consist of speckle 
observations obtained on UT 2010 June 18 with a two- 
color (approximately V and R) speckle c amera on the 
WIY N 3.5 m telescope on Kitt Peak (see iHowell et all 
|20H . and near-infrared ( J-band) adaptive optics (AO) 
observations conducted on UT 2009 September 8 with 
the PHARO camera on the 5 m Palomar telescope. No 
companions were detected around Kepler-10 within l'/5 
(for speckle) or 12'.'5 (AO), and more generally these ob- 
servations place strong limits on the presence of other 
stars as a function of angular separation (down to 0'/05 in 
the case of speckle) and relative brightness (companions 
as faint as AJ = 9. 5 for AO). These sen sitivity curves 
are shown in Fig. 9 of lBatalha et al.l (|2011| ) , and we make 
use of that information below. 

High-resolut i on spectra described also by 
IBatalha et ail (120111 ) and obtained with the HIRES 
instrument on the 10 m Keck I telescope place ad- 
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ditional limits on the presence of close companions 
falling within the spectrograph slit (Of! 87), such that 
stars within about 2 magnitudes of the target would 
generally have been seen. A small chance remains that 
these companions could escape detection if their radial 
velocity happens to be within a few kms -1 of that of 
the target (which is a narrow-li ned, slowly rotating star 
with vsini = 0.5 ± 0.5 kms" 1 ; IBatalha et all 1201 ID . so 
that the spectral lines are completely blended. This 
would be extremely unlikely for a chance alignment with 
a background/foreground star, but not necessarily for 
physically associated companions in wide orbits, i.e., 
with slow orbital motions. We explored this through 
Monte Carlo simulations. The results indicate that the 
probability of having a physical companion within a 
conservative range of ilOkms -1 of the RV of the target 
that would also go unnoticed in our speckle observations, 
and that additionally would not induce a RV drift on 
the target large enough to ha ve been detected in the 
high-precision measurements of [Batalha et al.l (| 2 1 If ) . is 
only about 0.1%. 

Finally, an analysis of the image centroids measured 
from the Kepler observations rules out background ob- 
jects of any brightness beyond about 2" of the target. 
This exclusion limit (equivalent to half a pixel) is con- 
siderably more cons ervative than the 0'.'6 reported by 
IBatalha et al.l ([201 ID , and accounts for saturation effects 
not considered earlier (given that at Kp = 10.96 the star 
is very bright by Kepler standards) as well as quarter- 
to-quarter variations (where "quarters" usually represent 
3-month observing blocks interrupted by spacecraft rolls 
required to maintain the proper illumination of the solar 
panels). 

4. APPLICATION OF BLENDER TO KOI-072.02 

The Kepler photo metry used her e is th e same as em- 
ployed in the work of IBatalha et al.l (|2011l ). and was col- 
lected between 2009 May 2 and 2010 January 9. These 
dates correspond to Kepler Quarter (first nine days of 
commissioning data) through the first month of Quar- 
ter 4. For this study we used only the long-cadence ob- 
servations (10,870 measurements) obtained by the space- 
craft at regular intervals of about 29.4 min. All blend 
models generated with BLENDER were integrated over this 
time interval for comparison with the measurements. 
The original data have been de-trended for this work 
by removing a first-order polynomial, and then apply- 
ing median filtering with a 2-day wide sliding window. 
Observations that occur during transits were masked 
and did not contribute to the median calculation. Be- 
cause this sliding window is considerably shorter than 
the 45.3-day orbital period, any ellipsoidal variations 
present in the original data should be largely preserved, 
although in any case they are expected to be very small 
for binaries with periods as long as this. We adopted 
also the e phemeris of mid-trans it for KOI-072.02 as re- 
ported bv IBatalha et al.l (|2011l) . which is T c [BJD] = 
2,454,971.6761 + Nx 45.29485 days, where TV is the num- 
ber of cycles from the reference epoch. 

Because it is relatively bright (Kp — 10.96), Kepler- 
10 was also observed by the Mission with a shorter ca- 
dence of approximately 1 min for a period of several 
months to allow an asteroseismic characterization of the 
star. A total of 19 oscillation frequencies were detected, 



and enabled a very precise determination of the mean 
stellar density. When combined with stellar evolution 
models and a spectroscopic determination of the effec- 
tive temperature and chemical composition, the result- 
ing parameters for the star are very well determined. 
Kepler-10 is relatively old (> 7.4 Gyr) but is otherwise 
quite similar to the Sun, with a temperature of T e ff = 
5627±44K, a mass and radius of M* = 0.895±0.060 M© 
and i?* = 1.056 ± 0.021 R Q , and a composition [Fe/H ] 
= -0.15±0.06 slightly below solar (jBatalha et al.ll20rl . 

As indicated earlier we considered four general sce- 
narios for false positives: chance alignments (a pair of 
background/foreground eclipsing objects) and hierarchi- 
cal triple systems, each with tertiaries that can be either 
stars or planets. The free parameters were varied over the 
following ranges: secondary mass Mi between 0.10 and 
1.40 Mq, in steps of 0.02 M©; tertiary mass M3 between 
0.10 and M2, also in steps of 0.02 M©; tertiary radius R3 
between 0.06 and 2.00i?j up in steps of 0.02i?j up ; impact 
parameter b between 0.00 and 1.00 in steps of 0.05; rela- 
tive duration D/D c i IC between 0.2 and 4.6 in steps of 0.2, 
corresponding to eccentricities up to 0.92 and values of 
to of 90° and 270° (see Sect.HJ); and relative distance A5 
(distance modulus difference) between —5.0 and +9.0 in 
steps of 0.5 mag, except for hierarchical triple configura- 
tions, for which A5 = 0. 

The goodness of the fit of each of the large number 
of synthetic light curves generated by BLENDER is quan- 
tified here by computing the \ 2 statistic and comparing 
it with that of the best planet model fit. The difference 
can be assigned a significance level (or false alarm rate) 
that depends on the number of free parameters of the 
problem. For example, for a blend scenario correspond- 
ing to a hierarchical triple system (4 degrees of freedom), 
a trial model giving a worse fit than the planet solution 
by A\ 2 = 4.72 is statistically diffe rent at the la le vel, 
assuming Gaussian errors (see, e.g. JPress et~aT1ll992l ). A 
fit that is worse by A\ 2 = 16.3 is different at the 3<r level. 
Hierarchical triple blends giving poorer fits than this are 
considered here to be ruled out by the Kepler photome- 
try. For background/foreground scenarios (5 degrees of 
freedom) the 3<r blend rejection level is A\ 2 = 18.2. 

4.1. BLENDER results 

In this section we describe the simulations carried 
out for the four general blend configurations mentioned 
above. Although the secondaries for the background sce- 
narios can in principle also be evolved stars (giants), as 
opposed to main-sequence stars, we consistently found 
that the transit light curves generated by such systems 
give a very poor match to the observations because they 
do not have the right shape (the ingress/egress phases 
are too long). Therefore, we restricted our exploration 
of parameter space to main-sequence stars only. 

An additional possibility for a false positive may stem 
from an error in the determination of the orbital period. 
If the true period were twice the nominal value, alternat- 
ing transit events would correspond to primary and sec- 
ondary eclipses, implicating a blended eclipsing binary. 
The primary and secondary eclipses would often (but not 
always) be of different depth. As part of the vetting pro- 
cess for each candidate, the Kepler Team examines the 
even-numbered and odd-numbered events to look for dif- 
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ferences in depth that may indicate a f a lse po sitive of this 
kind. As described by Bat alha et al.l ()201lD . no signifi- 
cant differences were found for KOI-072.02 beyond the 
2a level, where a represents the uncertainty in the tran- 
sit depth (9 ppm). Nevertheless, as the possibility still 
exists that the components of the eclipsing binary are 
identical, experiments were run with BLENDER to exam- 
ine the transit shape produced by such scenarios, and it 
was found that the ingress and egress phases are always 
much too long compared to the observations, as expected 
for two equal-size stars eclipsing each other. Thus, these 
scenarios are easily ruled out as well. 

4.1.1. Background eclipsing binaries (star + star) 

The simulations with BLENDER indicate that few back- 
ground blend scenarios with stellar tertiaries are able to 
mimic the transit features in the light curve at an ac- 
ceptable level, and they all correspond to somewhat ec- 
centric orbits. In Figure [5] we show the goodness of fit of 
these scenarios, with the small closed 3a contour repre- 
senting the region of parameter space within which the 
fits are satisfactory, according to the criteria given above. 
Only blends with secondary masses M% larger than about 
1.3 Mq are allowed, and the eclipsing binary can only 
be within a small range of distances behind the target 
(4.0 < AS < 4.7) for the dilution effect to be just right, 
such that the corresponding apparent brightness differ- 
ence AKp is between 2.5 and 3.5 mag (see figure). The 
best among these blend models (located near the center 
of the contour) provides a fit that is about 2.1c worse 
than a planet model (but still acceptable) , and is shown 
in the top panel of Figure [3] compared against the planet 
model. The tertiary stars in these blends are constrained 
to be very small, between 0.10 and 0.16 Mq. 

That most blends involving background eclipsing bi- 
naries can be ruled out may appear somewhat surpris- 
ing, and is worth investigating. Indeed, for a given 
measured transit depth dtran, a blend can only repro- 
duce the light curve if it contributes at least a frac- 
tion dtran of the total flux collected in the Kepler aper- 
ture. Thus, one would expect that binaries as faint as 
AKp = — 2.51og(d t ran) ~ 8.6 mag relative to the target 
should be able to match that amount of dimming if the y 
were totally eclipsed (see, e.g.. lMorton fc J ohnson 2011), 
and furthermore, that the measured duration could also 
be reproduced by a large range of secondary sizes with 
an appropriate combination of orbital eccentricity and u). 
Yet we find that no blends fainter than AKp = 3.5 give 
tolerable fits to the light curve (see Figure [5]). A visual 
understanding of the underlying reason for this may be 
seen in the bottom panel of Figure [3j in which we show 
a blend model that one would naively expect should be 
able to match the observations, according to the crude 
recipe described above. This particular blend scenario 
is marked with a cross in Figure [2J and corresponds to 
AS = 5 and M% — 1.0 Mq, resulting in a magnitude dif- 
ference of AKp = 5.6 for the EB relative to the target. 
While this model does yield a good match to the mea- 
sured depth, and even the total duration, it doesn't per- 
form nearly as well in the ingress/egress phases, which 
are too long when compared against the observations. 
The quality of this fit relative to the best planet fit, 
which can also be seen in the figure, corresponds to a 
lO.ler difference, and therefore BLENDER rejects it. Thus, 
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Fig. 2. — Map of the x 2 surface (goodness of fit) correspond- 
ing to a grid of blend models for KOI-072.02 involving background 
eclipsing binaries. The linear separation between the binary and 
the primary is cast in terms of the distance modulus difference. 
Contours are drawn as a function of the \ 2 difference from the 
best planet model fit (expressed in units of the significance level 
of the difference, a), and are plotted here as a function of the 
mass of the secondary star. Only blends within the small white 
contour yield acceptable fits to the light curve (within 3a" of the 
planet fit). Other colored areas correspond to regions of parame- 
ter space giving increasingly worse fits (4cr, 5cr, etc.), representing 
blends we consider to be ruled out. The Spitzer constraint is in- 
dicated by the shaded area: blend s with secondary masses in this 
region are excluded (see Sect. l3T2)l . although BLENDER itself already 
rules out all of these scenarios. Green lines running diagonally 
from the lower left to the upper right are labeled with the mag- 
nitude difference AKp of the blended binary relative to the target 
star. The hatched region below the AKp = 2 mag line represents 
blends with secondary stars bright enough that they would gener- 
ally be detected in our spectroscopy. Viable blends within the 3<r 
contour are seen to be confined to a narrow range of magnitude 
differences (2.5 < AKp < 3.5, dashed green lines). The dashed 
line at AKp = 8.6 indicates the envelope for the faintest blends 
that would be capable of reproducing the measured depth based 
on simple-minded estimates described in the text. As seen, BLENDER 
provides much tighter constraints than this. The cross corresponds 
to a blend model that gives the fit shown in the bottom panel of 
Figure [3] 

the reason most blends of this class can be ruled out is 
ultimately the high precision of the Kepler light curves, 
which provides a very strong constraint on the shape of 
the transit light curve, and in particular on the size ra- 
tio between the secondary and tertiary, which sets the 
duration of the ingress and egress phases. 

4.1.2. Background/foreground star + planet pairs 

There is a very broad range of blends consisting of 
a background or foreground star transited by a planet 
(as opposed to a star) that are found by BLENDER to 
give satisfactory fits to the data, as shown in Figure 21 
These viable blends occupy the area below the 3a con- 
tour represented with a thick white line. Secondary stars 
of all spectral types (masses) are permitted, in principle, 
although in practice other constraints described below 
eliminate a substantial fraction of them. All of these 
blends involve secondary+tertiary pairs that are within 
4 magnitudes of the target in the Kepler passband (di- 
agonal dashed line in the figure). The tertiary sizes in 
these blends range from 0.42i?j up to 1.84i?j up . 

Our Warm Spitzer observations set a lower limit of 
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Fig. 3. — Long-cadence Kepler observations for KOI-072.02 com- 
pared with two different blend models involving background eclips- 
ing binaries (red lines), and shown against the best fit planet model 
for reference (black line). All models are integrated over the du- 
ration of one cadence (29.4 min). Top: The blend model shown is 
the one giving the best fit for this type of scenario (2.1a difference 
compared to the planet fit). Bottom: Example of a blend model 
(indicated with a cross in Figure[2]| that illustrates the use of shape 
information by BLENDER in a case that would naively be expected 
to be a viable false positive scenario (see text). This particular 
scenario corresponds to a secondary mass of M2 = 1.0 Mq and a 
distance modulus difference of 5 mag relative to the target, giv- 
ing a brightness difference in the Kp band of 5.6 mag. Although 
it matches the depth and total duration of the transit, the ingress 
and egress phases are not well reproduced, so that the overall qual- 
ity of the fit is poor and the blend is ruled out at more than the 
10- a level. 

about 0.77 M© for the secondary masses of these blends, 
as described earlier; scenarios involving redder stars 
would result in transits at 4.5 /im significantly deeper 
than we observe (i.e., deeper than the measured depth 
+ 3cr). This exclusion region is indicated by the shaded 
area. Additionally, blends that are much brighter than 
AKp = 2 woul d most likely have b een detected spectro- 
scopically (see lBatalha et al.1l2011|) . so we consider those 
to be ruled out as well. We indicate this with the green 
hatched region in the lower right-hand side of the fig- 
ure. Finally, the colors of the background/foreground 
configurations simulated with BLENDER provide a further 
constraint which is represented by the blue hatched area 
on the lower left of the figure. This swath of param- 
eter space is excluded because the blends are signifi- 
cantly redder than the color index measured for Kepler- 
10 (Kp - K s = 1.465 ± 0.029), by more than three 
times the uncertainty in the observed index. As a re- 
sult of these complementary constraints, the only section 
of parameter space remaining for viable blends involving 
star+planet pairs is the area under the 3<t contour and 
limited from below and on the left by the hatched ar- 
eas (color and brightness conditions) and shaded area 
(Spitzer constraint), respectively. All of these blends 
have the eclipsing pair behind the target (foreground sce- 
narios are all ruled out). 

We note that in this star+planet blend scenario white 
dwarfs can also act as tertiaries, as long as they are 
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Fig. 4. — Similar to Figure [2] for blends involving background 
systems consisting of a star transited by a planet. The color scheme 
is the same as in Figure[2] Blends giving fits no worse than 3<r from 
the best planet fit are below the thick white contour labeled with 
that confidence level. The shaded left-hand side of the diagram 
corresponds to secondary masses excluded by constraints from the 
Spitzer observations. Lines of constant magnitude difference rela- 
tive to the target are shown in green, running diagonally from the 
lower left to the upper right. The dashed one at the top represents 
the boundary for the faintest viable blends (tangent to the white 
3ct contour). The solid green line below and parallel to it (AKp = 2 
mag) and hatched region to the right marks the area of parameter 
space excluded by our spectroscopic constraints. The blue curve 
and hatched region to the left represent blends that are excluded 
because they are too red in Kp — K s compared to the target. Note 
that the colors of the blended stars are computed from a different 
isochrone than that of the target, which explains why blends with 
secondaries of the same mass as the target are ruled out for being 
too red. The combination of all these constraints leaves only a re- 
duced area of parameter space (labeled "Allowed Region") where 
blend models give tolerable fits to the Kepler light curve and are 
not ruled out by any of our follow-up observations. 

cooler than the secondaries so that they do not lead to 
deep occultation events that would have been seen in 
the light curve of KOI-072.02. The above range of ter- 
tiary radii (0.42i?j up to 1.84i?j up ) excludes essentially 
all cool carbon-oxygen and oxygen-neon white dwarfs 
more massive than about 0.4 Mq, as these are smaller 
than the lower lim it set by BLENDER , which corresponds 
to 4.7 i?© (see, e.g. JPanei et al. 2000). Low-mass helium- 
core or oxygen-core white dwarfs that are the product of 
common-envelope evolution in binary stars can be con- 
siderably larger in size, although they appear to be very 
rare. The Kepler Mission itself has uncovered only thre e 
examples to date ()Rowe et al.ll2010HCarter et al.ll201lD . 
However, all of them are very hot (T e g > 10, 000 K), 
and produce deep and unmistakable fiat-bottomed oc- 
cupation si g nals. Model calculations such as those of 
iPanei et al.l (|2007| ) show that as these helium-core white 
dwarfs cool, their radii quickly become Earth-size or 
smaller. Therefore, we do not consider white dwarfs to 
be a significant source of blends for KOI-072.02. 



4.1.3. Hierarchical triple scenarios (star- 
star + planet blends) 



star and 



Eclipsing binaries composed of two stars physically as- 
sociated with the target are clearly ruled out by BLENDER, 
as they produce very poor fits to the Kepler light curves. 
For cases in which the tertiaries are planets, viable sce- 
narios identified by BLENDER span a range of secondary 
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Fig. 5. — Similar to Figure [2] for the case of hierarchical triple 
systems in which the secondary star is transited by a planet. The 
color scheme is the same as in Figure [2] In this case the vertical 
axis shows the tertiary sizes. The constraints from Spitzer (gray 
shaded area to the left of 0.77 Mq), color information (blue hatched 
area on the left), and spectroscopy (green hatched area on the 
right) are shown as in previous figures. When taken together these 
constraints eliminate all blends of this kind. 
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Fig. 6. — Same as Figure [5] but with the vertical axis showing 
the relative transit durations (D/D c ; rc ). 



masses and tertiary radii within the 3c contour shown 
in Figure [SJ Most of these configurations turn out to 
involve eccentric orbits, with transit durations longer 
than those corresponding to circular orbits along with 
secondary stars that are smaller than the primary (see 
Figure[6|). Once again other observational constraints are 
very complementary, and in this case they are sufficient 
to exclude all of these blends. For example, the shaded 
area of parameter space to the left of 0.77 Mq is elimi- 
nated by the Spitzer observations, as described earlier. 
The constraint on the Kp — K s color (hatched area on 
the left) is partly redundant with the NIR observations, 
but extends to slightly larger secondary masses. And fi- 
nally, the spectroscopic constraint removes the remaining 
scenarios corresponding to higher-mass (brighter) secon- 
daries. 



We conclude that of all the hierarchical triple blend 
scenarios that are capable of precisely reproducing the 
detailed shape of the Kepler transit light curve, none 
would have escaped detection by one or more of our 
follow-up efforts, including NIR Spitzer observations, 
high-resolution spectroscopy, or absolute photometry 
(colors). 22 This highlights the importance of these types 
of constraints for validating Kepler candidates, given 
that blends involving physically associated stars would 
generally be spatially unresolved by our high-resolution 
imaging with adaptive optics or speckle interferometry, 
and they would typically also be below the sensitivity 
limits of our centroid motion analysis, so that they would 
not be detected by those means. Therefore, the only 
blends we need to be concerned about for KOI-072.02 
are those consisting of stars in the background of the 
target that are orbited by other stars or by transiting 
planets. 

5. A PRIORI LIKELIHOOD OF REMAINING BLEND 
SCENARIOS FOR KOI- 72.02 

In order to estimate the frequency of the blend scenar- 
ios (i.e., background configurations) that remain possible 
after applying BLENDER and all other observational con- 
str aints, we follow a pr ocedure similar to that described 
bv iTorres etaU (|201lD for Kepler-9d. We appeal to 
the B esangon Galactic structure models of IRobin et al.l 
(2003) to predict the number density of background 
stars of each spectral type (mass) and brightness around 
Kepler-10, in half-magnitude bins, and we make use of 
estimates of the frequencies of transiting planets and of 
eclipsing binaries from recent studies by the Kepler Team 
to infer the number density of blends. Using constraints 
from our high-resolution ima ging (specifically, th e sen- 
sitivity curves presented by iBatalha et al.1 1201 ll their 
Fig. 9) we calculate the area around the target within 
which blends would go undetected, and with this the ex- 
pected number of blends. 23 

The recent release by iBorucki et al.l (|2011b[ ) of a list 
of 1235 candidate transiting planets (KOIs) from Kepler 
provides a means to estimate planet frequencies needed 
for our calcula tions, with significa nt advantages over the 
calculations of lTorres et alj (2011) for Kepler-9d, which 
were based on the ear lier list of candidates published by 
IBorucki et al.l (|2011a| ) . Not only is the sample now much 
larger, but the knowledge of the rate of false positives for 
Kepler is also much improved, and that rate is believed 
to be relatively small (20-40% dependin g on the level 
of vet ting of the candidate, according to IBorucki et al.l 
201 lat less than 10% acc ording to the recent study by 
Morton fc Johnson! [20ll . Thus, our results will not 
be significantly affected by the assumption that all of 

22 The possibility, remote as it may be, that the target has a 
physically associated companion that is nearly of the same bright- 
ness and that has managed to elude detection is always present 
(sec Sect. l3~3l) . not only here but in all previously discovered tran- 
siting planets. For the present purposes we do not consider this 
"twin star" scenar io as a false positive in the strict sense (see also 

ITorres et al. 2011), as the transiting object would still be a planet, 
only that it would be larger than we thought by about a factor of 
\/2 because of the extra dilution from the companion. 

23 In the case of KOI-072.02 the exclusion radius from the cen- 
troid motion analysis, 2" (Sect. [373) l. is significantly less constrain- 
ing than the high-resolution imaging, so is not as useful here as it 
was for Kcpler-9 d. 
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the candidates are planets (see also below). An addi- 
tional assumption we make is that this census is largely 
complete. Among these candidates we count a total of 
267 having radii in the range allowed by BLENDER for 
the tertiaries of viable blends (i.e., between 0.42 and 
1.84-Rjup)- W ith the total number of Kepler targets 
being 156,453 (jBorucki et al.ll2011bO . the relevant fre- 
quency of transiting planets for o ur blend calculation is 
/ planct = 267/156,453 = 0.0017. iSlawson et all (|20"ll 
have recently published a catalog of the 2165 eclipsing 
binaries found in the Kepler field, from the first four 
months of observation. Only the 1225 detached systems 
among these are considered here, since binaries in the 
category of semi-detached, over-contact, or ellipsoidal 
variables would not produce light curves with a shape 
consistent with a transit. The frequency of eclipsing bi- 
naries for our purposes is then /eb = 1225/156,453 = 
0.0078. 

Table [1] presents the results of our calculation of the 
frequency of blends, separately for background blends 
with stellar tertiaries (eclipsing binaries) and with plan- 
etary tertiaries. Columns 1 and 2 give the Kp magnitude 
range of each bin and the magnitude difference AKp rel- 
ative to the target, calculated at the upper edge of each 
bin. Column 3 reports the mean number density of stars 
per square degree obtained from the Besancon models, 
for stars in the mass range allowed by BLENDER as shown 
in Figure [2] In column 4 we list the maximum angu- 
lar separation p max at which stars in the corresponding 
magnitude bin would go undetected in our imaging ob- 
servations, ta k en fro m the information in the work of 
iBatalha et ail ([201 ll) . The product of the area implied 
by this radius and the stellar densities in the previous 
column give the number of stars in the appropriate mass 
range, listed in column 6 in units of 10~ 6 . Multiplying 
these figures by the frequency of eclipsing binaries /eb 
then gives the number of background star+star blends 
in column 7. A similar calculation for the background 
star+planet blends, making use of /planet) is presented in 
columns 7-10. We sum up the contributions from each 
magnitude bin at the bottom of columns 6 and 10. The 
total number of blends we expect a priori (blend fre- 
quency) is given in the last line of the table by adding 
these two values together, and is BF = 1.62 x 10~ 8 . 
The calculations show that background blends consisting 
of star+planet pairs contribute to this frequency about 
three times more than background eclipsing binaries. 

While we have assumed up to now that any com- 
panions to KOI-072.02 within AKp = 2 mag of the 
target would have been seen spectroscopically, we note 
that relaxing this condition to a much more conservative 
AKp = 1 has no effect at all on the contribution from 
eclipsing binaries, and a negligible effect on the contri- 
bution of star+planet scenarios. 

6. LIKELIHOOD OF THE PLANET INTERPRETATION FOR 

KOI-72.02 

To obtain a Bayesian estimate of the probability that 
KOI-072.02 is indeed a planet as opposed to a false pos- 
itive (or equivalently, the "false al arm rate", FAR) we 
follow the general methodology of iTorres et al.l ([201 ll ) 
and compare the a priori likelihoods of blends and of 
planets: FAR = BF/PF. If the a priori blend frequency 
is sufficiently small compared the planet frequency (PF), 



we consider the planet validated. Our a priori blend fre- 
quencies above correspond to false positive scenarios giv- 
ing fits to the light curve that are within 3c of the best 
planet fit. We use a similar criterion to estimate the 
a priori planet freque ncy by counting the KOIs in the 
(jBorucki et al.l 1201 lbl ) sample that have radii within 3cr 
of the best fit from a planet model (R p = 2.227±£o5? #e; 
see Tabled below). We find that 157 among the 1235 
KOIs are in this radius range (2. 06-2. 38 R®), giving PF 
= 157/156,453 = 0.0010. This results in a false alarm 
rate for KOI-072.02 of FAR = 1.6 x 10" 5 , which is so 
small that it allows us to validate the candidate with a 
very high level of confidence. The planet is designated 
Kepler- 10 c. 

This result rests heavily on the a priori frequency 
of planets from the Kepler Mission, derived from 
the assumption that all 1235 candidates reported by 
iBorucki et al.l (|2011b[) are indeed planets rather than 
false positives. If we were to be as pessimistic as to as- 
sume that as many as 90% of the small-size candidates 
are actually false positives (a similar rate of false pos- 
itives as is typically found in ground-based surveys for 
transiting planets), and at the same time that all of the 
larger-size candidates that come into the blend frequency 
calculation are real planets (thereby maximizing BF and 
minimizing PF), the false alarm rate would be 10 times 
larger than before, or 1.6 x 10~ 4 . This is still a very 
small number, and our conclusion regarding validation is 
unchanged. We note that a rate of false positives as high 
as 90% yields a planet frequency that is s trongly incon- 
sistent not only with the expectation s of IBorucki et al.l 
(|2011bfl and lMorton fc Johnsonl ([20111 ). but also with the 
independent results of ground ba sed Doppler surveys as 
reported by iHoward et aH ([2010D . 

In the above calculations we have implicitly assumed 
similar period distributions for planets of all sizes and 
for eclipsing binaries. However, it is conceivable that the 
results could change if the period distribution of plan- 
ets such as Kepler-lOc were significantly different from 
the one for larger planets that go into the blend fre- 
quency calculations, or from the one for EBs (which have 
a smaller contribution to BF; see Table [TJ. Therefore, as 
a further test we considered the impact of restricting the 
periods to be within an arbitrary factor of two of the 
Kepler-lOc period of 45.3 days, both in our blend fre- 
quency calculations and for the a priori estimate of the 
planet frequency, PF. We find that the planet frequencies 
are reduced by a factor of 4.5, and the eclipsing binary 
frequency by a factor of 10.4, and as a result the false 
alarm rate for KOI-072.02 is FAR = 1.4 x 10~ 5 , which is 
about the same as before. Thus, our conclusions are ro- 
bust against assumptions about the period distributions. 

Finally, our false alarm rate is conservative in the sense 
that we have not accounted for the flatness (coplanarity) 
of the Kepler-10 system. Only a small fraction of sin- 
gle transiting planets with periods as long as 45 days 
orbiting background stars (i.e., those acting as blends) 
are likely to transit, a priori, whereas a planet of this 
period such as Kepler-lOc is much more likely to tran- 
sit if it is coplanar with Kepler-10 b. Taking this into 
account would boost the planet frequency (PF) and de- 
crease th e FAR by as much as an order of magnitude 
(see, e.g. JBeattv &: Seag cr 2010). Coplanarity in multi- 
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pie systems is in fact supported by the large number of 
multiple transiting sy stem candidates found by Kepler 
(|Borucki et al.l2011btlLatfiam et al.ll2011l ) and their mu- 
tual in clinations seem to be small (1-5°: ILissauer et al.1 
12011 ri) . Therefore, we consider our estimate of the FAR 
for Kepler- 10 c to be conservative. 

7. DISCUSSION 

The stellar, orbital, and planeta ry parameters i n ferred 
for the system as determined by featalha eTaTl (|20ll 
are summarized in Table [2] to which we add the transit 
duration. The small formal uncertainty in the planetary 
radius (*~2.4%) derives from the relatively high preci- 
sion of the stellar radius, which is based on asteroseismic 
constraints on the mean density of the star. With its 
radius of about 2.2 i?®, Kepler- 10 c is among the small- 
est exoplanets discovered to date. The mass is undeter- 
mined as the Doppler signa t ure ha s not been detected. 
Nevertheless, iBatalha et al.l (|2011| ) placed a constraint 
on it based on the distribution of masses resulting from 
the Markov Chain Monte Carlo fitting procedure they 
applied to the existing radial-velocity measurements of 
Kepler-10. Their conservative 3-er upper limit for the 
mass is 20 M ffi . The corresponding maximum mean den- 
sity is 10 g cm^ 3 . 

Given a precise radius measurement and mass upper 
limit of 20 M s , some minimal constraints can be placed 
on the compo s ition of Kepler-10 c. Using the models of 
iFortnev et all (|2007l ). we find that an Earth-like rock- 
iron composition is only possible at ~ 20 M®. Lower 
masses would require a depletion in iron compared to 
rock, or more likely an enrichment in low-density volatiles 
such as water and/or H^/He gas. A 50/50 rock/water 
composition yields 2.23 R<§, at 7M®. Still lower masses 
are possible wit h a H^/He ga s envelo pe. Using mod- 
els presented in ILissauer et al.l (|2011al ) , a planet with a 
rock/iron core and a 5% H 2 /He atmosphere (by mass) 
matches the measured radius of Kepler- 1 c at only 3 Af© . 
A massive 20 M® core should have attained a H 2 /He en- 
velope, and it would appear to be stable at Kepler-10 c's 
relatively modest irradiation level, which would lead 
to a planetary radius dramatically larger than 2.23 i?©. 
This would tend to favor a s cenario where Kepler-10 c 
is more ak in to GJ 1214b (iCharbonneau et al.l 120091 : 
iBean et al.l 120101: iNettelmann et all I2010D and Kepler- 
11b and Kepler-llf, which are all below 7M® and en- 
riched in volatiles. 

The well measured inclinations of both Kepler-10 b and 
Kepler-10 c allow us to put a weak constraint on the true 
mutual inclination (</>bc) between the orbital planes of 
the two planets. Although the relative orientation in 
the plane of the sky (i.e., the mutual nodal angle) is 
unknown, the different impact parameters and result- 
ing apparent inclinations place a lower li mit on </>t, c . As 
discussed bv lRagozzine fc Holmanl (|2010t l. the geometric 
limits to the mutual inclination are given by \ib — i c \ < 
0bc < ib + i c , where i b = 84?4t^g (|Batalha et al.ll2011| ) 
and i c = 89?65±o ? 2 ) (Table [2]) are the usual inclinations 
with respect to the line of sight. Assuming a random ori- 
entation of the lines of nodes (which does not account for 
the a priori knowledge that both planets are transiting) , 
the mutual inclination is constrained to be in the interval 
5?25 < (f>bc < 174?05, with the most likely values being 



at the extremes of this distribution. Making the rea- 
sonable supposition of non-retrograde orbits, a mutual 
inclination close to the lower limit of about 5° is most 
likely for these planets. A more detailed probabilistic ar- 
gument requires making assumptions about the number 
of planets in the Kepler-10 system. 

This mutual inclination is on the high end of the distri- 
bution inferred for other Keple r multip le candidate sys- 
tems (1~5°) by ILissauer et al.l (|2011bl ). If this mutual 
inclination is typical for planets in this system, then it 
is relatively likely (depending on the orbital period) that 
other planets, if present, are not transiting. When con- 
sidering the set of Kepler candidates in multiple systems 
that have periods less than 125 days, the ratio of peri- 
ods between Kepler-lOc and Kepler-lOb (which is 54.1) 
is by far the highest of all period ratios of neighboring 
pairs of Kepler candidates (the next highest being 23.4), 
and is even higher than the period ratios between non- 
neighboring planets. Clearly, there is room for multi- 
ple additional planets between Kepler-10 b and Kepler- 
10 c. The preponderance of tightly-packed Kepler mul- 
tiple candidate systems suggests that additional planets 
may exist, and these may be revealed in the future with 
more detailed transit timing variation measurements. 

Kepler-10 c is the first Kepler target observed with 
Warm Spitzer with the aim of testing the wavelength 
dependence of the transit depth. This is currently the 
only facility available that has the capability of detecting 
such shallow transits at wavelengths that are sufficiently 
separated from the Kepler passband to be helpful. In 
this case the observations were successful, and the tran- 
sit at 4.5 /im is shown to have virtually the same depth 
as in the optical. This places a very strong constraint 
on the color of potential blends, which are restricted to 
have secondaries of similar spectral type as the primary 
star. 

The detailed analysis of the Kepler photometry with 
BLENDER combined with constraints from other observa- 
tions eliminates the vast majority of possible blend sce- 
narios. This includes most background eclipsing bina- 
ries (leaving only a small range of possible spectral types 
and relative fluxes for the secondaries), most of the sce- 
narios involving chance alignments with a star transited 
by a larger planet, and all possible hierarchical triple 
configurations. The latter are among the most difficult 
to detect observationally since they are typically spa- 
tially unresolved. The key factors that have allowed this, 
and made possible the validation of the planet, are the 
high-precision of the Kepler photometry, the relatively 
short ingress and egress phases (which places strong con- 
straints on the size ratio between the secondary and ter- 
tiary), and the near equatorial orientation, resulting in 
a relatively flat transit that leaves less freedom for the 
parameters of the eclipsing binaries. We expect BLENDER 
to be similarly effective for other Kepler candidates that 
show similar features in their light curves. 

Kepler-lOc along with Kepler-9d and Kepler-llg are 
examples of transiting planets that have not received the 
usual confirmation by dynamical means that previous 
discoveries have enjoyed (including essentially all ground- 
based discoveries), in which either the Doppler signature 
is detected unambiguously (and verified by the lack of 
bisector span variations), or transit timing variations in 
a multiple system are directly measured (as in Kepler-9 b 
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and c as well as the five inner planets of the Kepler- 11 
system). Instead, the planets in those three cases have 
been validated statistically, with a Bayesian approach to 
estimate the probability that the transit signals are due 
to a planet rather than a false positive. This probability 
has been computed by first estimating the a priori likeli- 
hood of a false positive, and then comparing it with the 
a priori chance of having observed a true planet. In the 
three cases mentioned above the ratio of the false positive 
to planet likelihoods is small enough that the planetary 
nature of the signal is established with a very high de- 
gree of confidence. For Kepler- 10 c the false alarm rate 
is 1.6 x 1(T 5 . 

The recent work of iMorton fc Johns on (201l]) has pro- 
vided a means of assessing a rough false alarm rate for 
Kepler candidates as a function of the depth of the tran- 
sit signal and the brightness of the object. As noted also 
by those authors, while these estimates are extremely 
valuable for statistical studies, the validation of candi- 
dates on an individual basis with a sufficiently high de- 
gree of confidence will usually require a much more de- 
tailed analysis of false positives, such as we have per- 
formed here. Masses for these objects (other than upper 
limits) may of course be difficult or impractical to deter- 
mine in many cases, but it is worth keeping in mind that 



some of the most exciting candidates to be discovered by 
Kepler will be in this category, namely, Earth-size plan- 
ets in the habitable zones of their parent stars. Except 
for stars of late spectral type, the RV signals will gener- 
ally be very challenging to detect with the sensitivity of 
current instrumentation. Thus, statistical validation of 
planets is likely to play an important role for Kepler in 
the years to come. 
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TABLE 1 

Blend frequency estimate for KOI-072.02. 
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Blend frequency (BF) = (0.0041 - 


h 0.0121) x 10" 6 = 


1.62 x 10 
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Note. — Magnitude bins with no entries correspond to brightness ranges in which all blends arc ruled out by a combination of BLENDER and other constraints. 
a The number densities in columns 3 and 7 differ because of the different secondary mass ranges permitted by BLENDER for the two kinds of blend scenarios, 
as shown in Figure \2\ and Figure HI 
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TABLE 2 

Star and planet parameters for Kepler-10 c. 



Parameter Value Notes 



Spectroscopically determined stellar parameters 

Effective temperature, T off (K) 5627 ±44 A 

Surface gravity, log g (cgs) 4.35 ± 0.06 A 

Metallicity, [Fe/H] -0.15 ±0.04 A 

Projected rotation, usini (kms - 1 ) 0.5 ± 0.5 A 

Inferred host star properties 

Mass, Mi, (M Q ) 0.895 ± 0.060 B 

Radius, R t (R Q ) 1.056 ±0.021 B 

Surface gravity, log (cgs) 4.341 ± 0.012 B 

Luminosity, L* (Lq) 1.004 ±0.059 B 

Absolute V magnitude, M v (mag) 4.746 ± 0.063 B 

Age (Gyr) 11.9 ±4.5 B 

Distance (pc) 173 ±27 B 

Transit and orbital parameters 

Orbital period, P (days) 42.294851 ^0076 C 

Mid-transit time, T c (HJD) 2,454,971.6761 Iq.oom C 

Scaled semimajor axis, a/i?* 49. 1 j^^'g C 

Scaled planet radius, R p /R* - 01938 t :oo024 C 

Impact parameter, b 0.299^q'q73 C 

Orbital inclination, t (deg) 89.65l|j'°2 C 

Transit duration, A (hours) 6.863^Q Qgg C 

Parameters for Kepler-10 c 

Radius, R p (R @ ) 2.227±°' u ^ B,C 

Mass, M p (M©) < 20 D 

Mean density, p p (g cm -3 ) < 10 D 

Orbital semimajor axis, a (AU) 0.2407t°'oos3 B 

Equilibrium temperature, T eq (K) 485 F 



Note. — In most cases these parameters are taken from Batalha ct al. 
(2011). A: Based on an an alysis by D. Fischer of the Keck/HIRES templat e 
spectrum using SME (sec Valcnti & Piskunov 1996; Batalha ct al. 2011); 
B: Based on the astcroscismology analysis and stellar models; C: Based on 
an analysis of the photometry; D: Upper limit corresponding to three times 
the 68.3% credible interval from the MCMC mass distribution; E: Based 
on Newton's revised version of Kepler's Third Law and the results from D; 
F: Calculated assuming a Bond albedo of 0.1 and complete redistribution 
of heat for rc-radiation. 



